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5 Prinzipien der Rontgen-Cm 



■tomographic 



5.1 Prinzipicn 

5.1.1 Weichteilkontrast bzw. Stralilcnschuacliung 
bei dcr klassischcn Rontgcntcchnik 

Oblichc Rontgenaufnahmen bilden die zwctdimcnsionale Projektion eines drei- 
dimensionalen Korpers ab. wobci jcder Bildpunktwert dicser Abbildung die 
Schwiichung seines zugeordneten Strahls Kings des Strahlcnweges durch den 
Korper wiedergibt. Kleine Unterschiedc im Slrahlenschwachungsvermogen 
etwa von Weichteilgewcbe innerhalb des Korpers kommen gewohnlich kaum, 
bei gleicher Schwiichung benachbartcr Strahlbcrciche iiberhaupt nicht zur Gel- 
tung. Wenn z. B. ein Korpersegment dcr Dicke d x mil einem Schwachungskoeffi- 
zienten neben einem zweiten der Dicke d 2 und dem Koeffizienten pi 2 
durchstrahlt wird, und wenn ^ l il l =fi 2 'd 2 sein sollte. wird kein Unterschied 
zu beobachten sein. 

Alfgemein wird die Schwiichung durch Linienintcgrale 

s 

j 

iiber die Verteilung der Schwiichungskoeffizienten //(.v. y, i) des betrachteten 

MeBobjekts liings der Strahlenwege a? bestimmt. Bei einem realen Korper ist 

es nicht unwahrscheinlich, daB benachbarte Linieninregrale sich kaum unter- j 

scheiden. 

In der konventionellen Rontgendiagnostik behilft man sich vielfach mit der 
Applikation unterschiedlicher Kontrastmittei, urn das Strahlenschwachungsver- 
mogen bestimmter Bereiche anzuheben, soweit dies klinisch angezeigt ist. Unver- 
meidlich bleibt, daB wegen der Uberlagerung aller Strukturen liings der Strah- 
lenwege interessierende Areale nicht oder nicht eindeutig lokalisiert werden 
konnen. Auch die klassische Verwischungstomographie setzt die Lokalisierung 
eher voraus als sie zu liefern. \ 



(5.1) 



5.1.2 Einfuhrung in das Verfahren 

Die Rontgen-Computertomographie - kurz mit CT bezeichnet - ist ein speziel- 
les Rontgen-Schichtaufnahmeverfahren, das sich im Bijdaufbau grundsiitzlich 
von dem klassischen, nach dem Verwischungsprinzip arbeitenden Rontgen- 
Schichtaufnahmeverfahren [5.3] unterscheidet. Bei CT-Aufnahmen ergeben sich 
Transversalschnittbilder, also Abbildungen von Korperschichten. die im wesent- 
lichen senkrecht zur Korperachse orientiert sind. Diese neue Schichtdarstel- 
lungsmethode stellt nun - im Gegensatz zur klassischen Rontgentechnik - die 
Verteilung des Schwachungswertes p(x, y, z) selbst dar. Das gibt ihr das Poten- 
tial, interessierende Organe prinzipiell in einem dreidimensionalen Bildraum 
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.5./ i'rinzipirn 
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darzustellen und geringe Dichtcuntcrschicde von Weichieilgewebc gut erkenn- 
bar. d.h. kontrastreich abzubifden. 

Allerdings benotigt die CT zur Losung dieser Aufgabe Projeklionen unter schr 
vicfen Winkeln. also ungieich vie! mehr information, als bci der dircktcn Projek- 
lionsabbildung der kiassischen Rontgentechnik genutzt wird. Sic reduzieri diese 
Informationsmenge aber auch gleich wieder, indem sic sich primar auf die Dar- 
stcllung eines nur zweidimensionalen Schnittes durch den Korper beschrankt: 
Ansteilc von //(.v, y, z) wird nur // z (.v, y) in der Position z (Position dor Schicht 
auf der Korperlangsachse) bestimmt. Aus mef3technischen Griindcn wird dieser 
Schnitt auf cine schmaie Schicht geweitet, deren Dicke h gewohnlich wahlbar 
gehalten, und uber die wiihrend der Messung gemitteh wird. Schichtdicken 
von 1 mm bis 10 mm sind derzeit gebrauchlich und wiinschenswert. Die Schnitte. 
als deren geometrischer Ort die Mittellinien der Schichten anzuschen sind, liegen 
ublicherwcise in einem Winkelbereich von ±30° um die Scnkrcchtc zur Langs- 
achse des Paiionten, weshalb man auch zuweilen von „transaxialer CT" sprichl. 

Im Begrifl* Tomographic kommt der Schnittbildcharakter des Darstellungsver- 
fahrens zum Ausdruck, das Attribut Computer deutet auf das rechnende Ele- 
ment hin, das zur Bestimmung von /i(.v, y) zwangsliiufig benotigt wird. Die 
vollstiindige dreidimensionale Rekonslruktion von /i(.v, y. z) gelingt mittelbar. 
indem man nacheinander Schicht an Schicht setzt und venniBt. Dieses schritt- 
weise Vcrfahren wird jedoch seit 1989 mehr und mehr durch die sogenannte 
„Spiral-CT" abgelost, bei der die Verschiebung des Patienten in Korperlangs- 
richmng kontinuierlich wahrend der MeBdatenerfassung geschieht (s. Ab- 
schnitt 5.2). 



5.1-3 McIJprinzipien 

■ 

Es ist naheliegend. zur Erzeugung einer Schichtaufnahme den von der Rontgen- 
rohre emittierten Strahlenkege) so auszubJenden, daB ein ebener StrahJenfacher 
entsteht, der eindimensionale Zentralprojektionen der durchstrahlten Schicht 
entwirft. Zur exakten Rekonstruktion der Verteilung der Schwachungswerte 
// z (.\\ v) muB dieser Strahlenfacher senkrecht auf der Drehachse stehen und 
auBerdem so weit gespreizt sein, daB er aus jeder Projektionsrichtung die anvi- 
sierte Schicht des MeBobjekts vollstandig iiberdeckt. Computertomographen. 
die soiche Strahlenfacher nutzen, sind die derzeit gebrauchlichsten CT^Gerate. 
Es gibt sie in zwei Varianten: 

O In einer ersten Variante. bei den sogenannten „FacherstrahlgerateiV (auch: 
..direkte; Facherstrahlgerate' 1 ), rotieren Rohre und Strahlenempfanger 
gemeinsam um eine Drehmitte, die auch Mitte des kreisformigen MeBfeldes 
ist(Bild 5.2c). 

t> In der zweiten Variante. bei den sogenannten „RingdetektorgeraterT (auch: 
..inverse Facherstrahlgerate"), rotiert die Rohre allein innerhalb eines Rings 
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5 Prinzip'wn der Ri'mt^cM-CompuWrton 




o - o - 




c) " ' d) 



a) Translation-Rotations-Scanner mil cinem Strahlcnbundcl und einem Dctcktor 

b) Translation-Rotations-Scanner mit mchrercn Strahlenbiindeln und mchrcrcn Dctckto- 
ren (hicr nur zwei gezeichnet) 

c) Facherstrahlgcrai mit rotierendem Dcteklorsystcm (direktes Facherstrahlgcrai) 

d) Ringdetckiorgeral mit stalionarem Detektorsystem (inverses Facherstrahlgcrai) 

Bild 5.2 Unterschiedlichc Abtastprinzipien von Compulcrlomographcn 



von Detektorelementen, von denen jedes als Konvergenzpunkt eines sich 
wahrend der Rotation entwickelnden Strahlenfachers betrachtet werden 
kann (Bild 5.2d). 

Bei der ersten Variante bestimmt die aktuelle, mittlere Fokuslage, unter der 
MeGvverte aus dem Strahlenempfanger ausgelcsen werden, den Winkel einer 
Projektion. Bei der zweiten Variante legt die Position eines jeden Detektorele- 
mentes auf dem Detektorring den jeweiligen Projektionswinkel Test. 

5.1.4 Bilderzeugung 

Das Bildrekonstruktionsverfahren laBt sich am besten an einem. MeBprinzip 
erlautern, das Parallelprojektionen liefert und in den .,Translations-Rotations- 
Scannern" realisiert wurde (Bilder 5.2a und b). Diese auch Linearscanner 
genannten CT-Gerate standen am Anfang der CT-Entwicklung. 

Wie in Bild 5.3 dargestellt, wird hier ein einzelner Strahl der Breite b und Hohe 

(Dicke) h mit einem Blendensystem (Koilimator) von der Rontgenrohre durch ; 

das MeBobjekt hindurch zu einem Detektor gefuhrt. Mit der Hohe h ist die" 

Schichtdicke festgelegt, die Breite b begrenzt das raumliche Auflosungsvermogen j 
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Bild 5.3 

MeBanordnung des einfachste 
(Linearscanner nach Bild 5.2 a 
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5.1 Prinzipien 



der McBcinrichtung. Rohrc und Dctcktor sind uuf einem gcnicinsumen Tragcr 
(„Gantry") monticrt, lassen sich senkrecht zum Strahl vcrschicbcn und um die 
MeBfeldmitte drchen. Der Detcktor reagiert auf die eintreffenden Strahien mit 
elektrtschen Signalen /, deren' Amplitude proportional zur Intensilat dicser 
Strahien ist. Aus dem Verhaltnis dieser Signale zum crwartetcn Signal / 0 ohne 
schwachendcs Objekt ergeben sich die gerade vorliegenden Strahlcnschwachun- 
gen durch das MeBobjekl. 

Auf die durchstrahlte Schicht sei ein ortsfestes Koordinatcnsystem (x, y) gclegt, 
in dem die Objektfunktion /i(x, y) lokalisiert ist, ferner ein zweites, zum ersten 
konzentrischcs System (c, y //), dessen c-Achse parallel zum MeUstrahl ausgerich- 
tet ist, das also die Drehung der Gantry mit vollzicht (Bild 5.4). In jenem (f, #/)- 
System registriert der Detektor die Intensitatsprofile /(</>, j/), wenn Rohre und 
Detektor - genauer, wenn der Strahl - parallel zu sich selbst verschoben, nach 
erfolgter Traverse um einen klcinen Winkel A(p gedreht, dann riicklaufig ver- 
schoben und wieder gedreht wird. Dieses Wechselspiel erfolgt so langc, bis . 
ein aufgelaufener Drehwinkel von 180° erreicht ist. Zu jedcm Winkel <f> gibt 
es danach einen Satz paralleler Strahien, deren Signale in /(</>, ;/) gesammelt 
wurden. 

Die Schwachungsprofile 

. *)// 0 (<M) (5.2) 

biiden also Parallelprojektionen der Schicht zum Projektionswinkel (p. Da es 
gleichgultig ist, von welcher Seite ein Objekt durchstrahlt wird, sind Schwii- 




Bild 5.3 

Mefianordnung des einfachsten Translations-Rotations-Scanners 
(Linearscanner nach Bild 5.2 a) 
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5 I'rinzipien tier liihUgc.n-Coniputcn^jBgrtiphic 




p(4>, ;?) = Jm(|. 



Qetektor 



x. y Raumfeste Koordinatcn 
<;, // An das Mefosyslem 

gcbundenc Koordinatcn 
(f) Projektionswinkei 



/ Ronlgcnsirahlinlcnsiiat 
p SchwachungskocfTizicni 
s Weg cincs Slruhls 



Bild 5.4 Winkcl- und Ortskoordinatcn dor Parallcl.strahlgeomctrie 



chungsprofile. die unter den Winkeln <f> gemesscn werden, identisch mit dencn, 
die bei 04-180° zu gewinnen sind; bei Gleichsetzung ist nur aufdie gegenliiufige 
Ausrichtung von jj zu achten, d. h., es ist 



S(ct> + n. n) = S{<t>. -17) 



(5.3) 



anzusetzen. 

Wenn man noch annimmt, die Strahlen batten einheitlich erne einzige Energie 
(z.B. 70keV), so da(3 alien Substanzen, mit denen die Strahlung in Wechselwir- 
kung tritt, ein bestimmtes spezifisches Strahlenschwachungsvermogen /* (d.h. 
Schwachungskoeffizient) zugeordnet werden kann, dann schwiicht sich die 
Anfangsintensitat I 0 {(p. >/) ab auf (s. Abschnitt 2.1.2) 



x = C cos </> — /? sin<p 
y = c sin <p + '/ cos (j> 



(5.4) 



Bildet man die Logaiithmen, so erhalt man die Linienintegrale der Projektionen 
p(^/;) = Jn/ 0 /7 = J/i(.Y,>')dc. (5.5)' 



deren Gesamtheit auch als die Radon-Transformierte der Verteilung /i(x, y) 
bezeichnet wird (s. Abschnitt 2.3.1). 



Wie in Kapiiel 2 erl 
aus /?(</>, }j) crrechenl 
(s. AbschniU 5.4.1), b 
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Wie in Kapilcl 2 erliiutert, ist dicse Transformation umkchrbar, fx{x, y) also 
aus p{(/>, i]) crrcchcnbar. Gebriiuchlich ist der sogenannto Fallungsalgorithmus 
(s. Abschnitt 5.4.1), bei dcm die Linienintegrale jc Projektion zunach'st mit einer 
spcziellen Funktion („Faltungskern" genannt) gefaltet und dann langs der 
urspriinglichcn Strahlrichtungen auf die Bildebcne riickprojiziert werden 
(Bild 5.5). Dieses Verfahren hat den Vorzug, daB die notwendigcn Rechenpro- 
zesse im Prinzip schritthaltend zur Messung durchlaufen werden konnen - aller- 
dings unter der Voraussetzung, daB hinreichend Icistungsfahigc Prozessoren 
zur Verfiigung stehen. Dies war z. B. kein Problem bei den crstcn Linearscan- 
nern, deren Mef3zeit etwa 5 min je Doppelschicht betrug (gleichzcitige Messung 
mit zwei urn die Schichtdicke h versetzten Detektoren). Erst bei den Nachfolgc- 
geriitcn, bei denen mit MeBzeiten von etwa 20 s je Schicht oder auch 60 s je 
Doppelschicht mchrere Detektoren mehrere Parallelprojcktioncn je Traverse 
zuglcich aufnahmen (Bcispiel SIRETOM2000. Prinzip s. Bild 5.2b), muBien 
damals schon spezielle Prozessoren entwickclt werden. . 

Ein ganz ahnlichcs Verfahren ist auch bei den Facherstrahlgeraten an wend bar. 
Es gclang namlich, fur Projektionen, die nach Projeklionswinkel <p (d.h. nach 
dem Azimut des Fokus) und nach Strahlenwinkel ft innerhalb des Fiichers 
parametrisiert sind - also p(<p, fi) anstelle von /){</>, ;/) - ebenfalls einen analogen 
Faltungsalgorithmus anzugebcn (s. Abschn. 5.4.2), der auf die gleichc Prozessor- 




Projektion eines homogcnen Zyiinders 





Riickprojektion ohne Faltung 



Riickprojektion nach Faltung 



Bild 5.5 

Bilderzeugung durch Riickprojektion und Bedeutung der Faltung der MeGwertprofiie 
mit dem Faltungskern 
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5 Prinzipicn tier Rontt'i'ij-Cuni^fjFomographie 



struktur abbildbar ist (s. Bild 9.21). Allcrdings crzwingt hicr die Riickprojektion 
eincn wcsontlich hdhcrcn Rechcnaufwand, weshalb audi Transformationen von 
Facher- auf Parallelprojektionen vorgenommcn wcrden, wenn die streng schritt- 
haltcnde Berechnung nicht gefordert ist. Dies erfolgt entweder schon bei der 
Datengcwinnung odcr sputer durch Umsorticrung der Dalcn und spezielle Inter- 
polationen. Solche Umordnungen sind audi dann erforderlich, wenn Rekon- 
struktionsalgorithmen angewandt wcrden sollen, die zwangslaufig auf Parallel- 
daten aufsetzen -erwahnt sei nurdas direkte Fouricr-Rekonslruklionsverfahren 
(s. Abschnitt 5.4.1), das wcgen der Vcrfiigbarkcit spczicllcr Prozessoren eine 
gewisse Attraktivitat genicBt. 



5.2 Realisierung I 

5.2.1 Grcnzen des Strahler-Mettsystems | 

Von der grundsatzlichen Idee bis zur Verwirklichung praktisch einsetzbarer j 
CT-Gerate gibt es verschiedene Wege, die aber alle zu einigen unvermeidbaren j 
Einschrankungen zwingen. j 

i 



Abtastdichte ■ 

Eine erste Einschrankung betrifft die Abtastdichte: Weder bei den Parallelstrahl- 

geraten noch bei den Facherstrahlgeraten konnen MeBwerte in beliebig dichter j 

Folge abgelesen werden. Die in einer vorgegebenen Zeit maximal aufnehmbare ! 

Anzahl der MeBwerte je Projektion ist naturgemaft begrenzt. Dies gilt auch : 

fur die Anzahl der Projektionen je Umlauf (oder ,.Scan k \ wie der Vorgang der \ 

MeBwertaufnahme auch genannt wird). Ob sich diese Begrenzung ungunstig 

auswirkt, hangt wesentlich vom raumlichen Auflosungsvermogen der Anlage 

ab. In der Regel werden umso mehr MeBwerte benotigt, je hoher das Auflo- ! 

sungsvermogen ist. • i 

Auflosung 

Eine zweite Einschrankung betrifft das Auflosungsvermogen: Sowohl der Fokus 
der Rontgenrdhre als auch die Detektorelemente der Strahlenempfanger sind 
von endlicher GroBe. Der Fokus, weil er hohe thermische Energie aufnehmen 
muB, um eine geniigend groBe Strahlendosis emittieren zu konnen, die Detek- 
torelemente, weil sie einen hinreichend groBen Bruchteil der auftrefFenden | 
Strahlenmenge absorbieren mussen, damit brauchbar hohe Signale entstehen. f 
Die Verbindungslinien von Fokus und Detektorelementen definieren ein Strah- ! 
lenbiindel endlicher Breite, unter der das MeBobjekt abgetastet wird. I 



Detcktoren bei Fat 

Bei direkten Fiich 
Zahl der Detektoi 
350 und 800. Diesc 
daB die MeBfelder 
von etwa 45 bis 5( 
daraus, daB der F; 
nur zur Halfte be 
mit einem Halbfac 
auf der eihen Seite 
der Drehung jewei 
man Teilfacher ge 
ist beispielsweise e 

Ein weiterer Grui 
die das Gerat eing^ 
gen verlangt ist. b 
so betont wird, set. 



Detektoren bei Rin 

Bei den Ringdetek 
durch den zeitliche 
sind hier zwischer. 
nicht bekannt sin. 
die Anzahl der Pre 
gerat, wo die Anza 
big festgelegt were 
beschrankt wird. 

In diesen Gegebe; 
Anzahl der Detek 
sowohl durch Ani 
durch moglichst h 
durch die Belange 
eine stattliche Anz. 
ring. Das sind me) 
rechnung eigentlic: 

Dieser Aufwand li 
tionen uber 360° \ 
nen aus einem ^ 
Bedeckung des Mv 
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5.1,4 Image generation 

The image reconstruction process can best be explained by way 
of a measuring principle which provides parallel projections 
and was implemented in the "translation-rotation scanners" 
(Figs. 5.2a and b) - These CT devices which are also referred 
to as linear scanners marked the beginning of the CT develop- 
ment . 

As is illustrated in Fig. 5.3, a single beam of width b and 
height (thickness) h is guided by means of a diaphragm system 
(collimator) from the X-ray tube through the measuring object 
to a detector- Height h establishes the layer thickness, width 
jb defines the three-dimensional resolution capability of the 
measuring equipment * The tube and the detector are mounted on 
a common carrier ("gantry"), can be moved perpendicular to the 
beam, and can be rotated about the measuring field centre. The 
detector responds to the incoming beams with electrical sig- 
nals I, the amplitude of which is proportional to the intens- 
ity of these beams. The currently existing beam attenuations 
due to the measuring object are obtained from the ratio of 
these signals to the expected signal I 0 without attenuating 
object. 

A stationary coordinate system (x, y) be placed on the X-rayed 
layer, in which the object function p(x, y) is localized; fur- 
thermore, a second system r\) which is concentric to the 
first one, the | axis of which is oriented parallel to the 
measuring beam, so that it follows the rotation of the gantry 
(Fig. 5.4), In this q) system the detector records the 
intensity profiles tj) when the tube and the detector - 

more precisely, if the beam - is moved parallel to itself, is 
rotated through a small angle Atf after carrying out a traverse 
movement, then moved in the opposite direction, and then ro- 
tated again. This cyclic sequence is repeated until a summed 
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up angle of rotation of 180 Q is obtained. For each angle $ 
there is a set of parallel beams, the signals of which were 
collected in I{<p r rj) , — 




Measuring arrangement of the simplest translation-rotation 
scanner (linear scanner according to Fig* 5.2a) 



Thus, the attenuation profiles 

r\) = f?)/l 0 ($, tfl (5.2) 

foria parallel projections of the la) " to the projection angle 
0. Due to the fact that it is irrelevaj. . '* from which side an 
object is X-rayed, attenuation profile? . xich are measured 
under the angles 0 are identical with those which can be ob- 
tained at $ + 180°; upon equating, though, the opposite 
orientation of 7 must be taken into consideration, i.e. 



(5.3) 
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is to be written. 



x ' y 



i 



stationary coordinates 

coordinates associated with the measuring system 
projection angle 
X-ray intensity 
attenuation coefficient 



s 



travel of a beam 



Fig. 5,4 Angle and position coordinates of the parallel beam 



If one further assumes that the beams homogenous ly have a 
single energy (e.g. 70 kev) so that all substances with which 
the radiation interacts can be allocated a certain specific 
beam attenuation capability p (i.e. attenuation coefficient), 
then the initial intensity i 0 ($, ij) decreases to (see para. 
2.1.2) 



When forming the logarithms, one obtains the line integrals of 
the projections 



geometry 



(5.4) 



x = 4 cosf - t] sin0 
y ~ £ sin$J - tf co50 




(5.5) 



which, in their entirety, are also referred to as the Radon 
transform of the distribution jj(x, y) (see para. 2.3.1). 
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As has been explained in Chapter 2, this transformation is 



the so-called convolution algorithm is conventionally used 
(see para. 5.4.1)/ where the line integrals per projection are 
first convoluted by a special function {known as "convolution 
core") and then re-projected onto the image plane along the 
original beam directions (Fig. 5.5). This method is advantage- 
ous in that th^ necessary calculation processes can princip- 
ally be executed in step with the measurement - under the pre- 
condition, however, that sufficiently powerful processors are 
available. This did not present a problem e»g» with the first 
linear scanners whose measuring time amounted to approx. 
5 min. ptir double layer (simultaneous measurement by means of 
two detectors offset by the layer thickness h) . It was only 
with the successor devices, where several w tectors recorded 
several parallel projections per traverse simultaneously with 
measuring times of approx. 20 sec per layer or 60 sec per 
double layer (example SIRETOM 2000; principle see Fig. 5.2b), 
that already then special processors had to be developed. 

A very similar method can also be employed with the fanned 
beam devices* One succeeded in specifying an analogous convol- 
ution algorithm, too, (see para. 5.4.2) for projections which 
are parameterised to projection angle <p (i.e. to the azimuth 
of the focus) and to beam angle within the fan - i.e. p{<p,fi) 
instead of p($, ij) , which can be mapped on the same processor 
structure (see Fig. 9.21). 



reversible, i.e. p(x, y) can be calculated from 
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Projection of a homogenous cylinder 

Reverse projection without convolution 

Reverse projection after convolution 

Fig. 5.5 

Image generation by reverse projection and significance of the 
convolution of the measuring value profiles with the convolu- 
tion core 



The reverse projection, however, makes considerably higher 
calculation efforts mandatory, which is the reason why trans- 
formations from fanned to parallel projections ,-ire made if the 
strict in-step calculation is not required. This is done 
either as early as during data acquisition or later by data 
re-sorting and special interpolations. Such re-arrangements 
are also required where reconstruction algorithms are to be 
employed which inevitably are based on parallel data - with 
reference being made here to the direct Fourier reconstruction 
method (see para. 5-4.1) which, due to the u ' " " Mlity of 
special processors, enjoys a certain attractive le ,s. 
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